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ABSTRACT 

We report on the first Cosmic Origins Spectrograph (COS) observations of damped and sub-damped 
Lyman-a (DLA) systems discovered in a new survey of the gaseous halos of low-redshift galaxies. From 
observations of 37 sightlines, we have discovered three DLAs and four sub-DLAs. We measure the 
neutral gas density fin i, and redshift density cLV /dz, of DLA and sub-DLA systems at z < 0.35. 
We find dA/"/dz=0.25tg i4 and ft H i =l-4±^;f x 1CT 3 for DLAs, and cW/dz=0.08t^ with ft H i 
=4.2^g'g x 10~ 5 for sub-DLAs over a redshift path Az = 11.9. To demonstrate the scientific potential 
of such systems, we present a detailed analysis of the DLA at z Q b s =0.1140 in the spectrum of SDSS 
J1009+0713. Profile fits to the absorption lines determine log N(H I)=20.68 ± 0.10 with a metallicity 
determined from the undepleted element Sulfur of [S/H]=— 0.62 ± 0.18. The abundance pattern of 
this DLA is similar to that of higher z DLAs, showing mild depiction of the refractory elements Fe 
and Ti with [S/Fe]=+0.24±0.22 and [S/Ti]=+0.28±0.15. Nitrogen is underabundant in this system 
with [N/H]=— 1.40 ± 0.14, placing this DLA below the plateau of the [N/a] measurements in the 
local Universe at similar metallicities. This DLA has a simple kinematic structure with only two 
components required to fit the profiles and a kinematic width of Awgo=52 km s . Imaging of the 
QSO field with HST/WFC3 reveals a spiral galaxy at very small impact parameter to the QSO and 
several galaxies within 10", or 20 comoving kpc at the redshift of the DLA. Followup spectra with 
LRIS on the Keck telescope reveal that none of the nearby galaxies are at the redshift of the DLA. The 
spiral galaxy is identified as the host galaxy of the QSO based on the near perfect alignment of the 
nucleus and disk of the galaxy as well as spectra of an H II region showing emission lines at the QSO 
redshift. A small feature appears 0.70" from the nucleus of the QSO after PSF subtraction, providing 
another candidate for the host galaxy of the DLA system. Even with these supporting data, we are 
unable to unambiguously identify the host galaxy of the DLA, exemplifying some of the difficulties in 
determining DLA hosts even at low redshift. 
Subject headings: none yet 



1. INTRODUCTION 

Quasar (QSO) absorption line systems provide a 
unique means to study the intergalactic medium, as well 
as the interstellar medium (ISM) of galaxies at all red- 
shifts. Resonance absorption lines from metallic ions in 
the ISM of distant galaxies along the line of sight to a 
QSO provide a powerful method for the determination 
of chemical abundances at high redshifts, and also reveal 
the properties of low-density gas that is difficult or even 
impossible to detect with any other method, even in the 
nearby Universe. As this method is independent of the 
luminosity of the intervening galaxy, it is a relatively un- 
biased tracer of the chemical evolution of the Universe 
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as a whole. 

The H I A 1215 absorption line in Damped Lyman- 
a systems (DLAs, log N(H I)>20.3) and sub-DLAs 
(19.0<log N(H I)<20.3) can be fit for accurate col- 
umn density measurements due to the extended damp- 
ing wings in the profile. DLAs and sub-DLA sys- 
tems have been shown to contai n the major i ty of 
the neutral gas in the Universe dWolfe et al.1 119951: 
Peroux et all |2003bt [Prochaska. Herbert-Fort fe Wolfd 
2005HNoterdaeme et al.ll2009D . 

The DLAs are of particular importance for galactic 
chemical evolution studies. With such high column den- 
sities, the absorption lines from metal atoms are easily 
visible. Also, systems with high column densities are ex- 
pected to remain mainly neutral due to self shielding, 
alleviating the need for uncertain ionization corrections 
to the abundances. 

There is still much debate as to the nature of the galax- 
ies hosting DLAs and sub-DLAs. Even at low redshift, 
finding a faint galaxy near the bright point source of the 
QSO is challenging. At higher redshift, the cosmologi- 
cal dimming of surface brightness (fj, oc (1 + z) A ) makes 
galaxy detections even more difficult. Nonetheless, these 
systems are a unique laboratory to study the ISM of 
galaxies over a wide range of redshifts. 

At z > 1.65, the Lyman-a line shifts above the atmo- 
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J0925+4004 log N(H I)= 19.55±0. 15 z ab! =0.247 



10928 + 6025 log N(H I)= 19.35±0. 15 2,^=0.1538 
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Figure 1. Lyman-a absorption lines from the DLAs and sub-DLAs in this sample. The solid green line is the adopted continuum, and 
the best fit profiles are bracketed by profiles with column densities la above and below the best fit value. The lcr flux uncertainty is shown 
in cyan below the data. 



spheric cutoff at ~ 3000A, and ground based surveys 
can efficiently search for DLA systems. Severai such 
surveys have determ i ned th e z > 2 neutral gas density 
ftHi (IWolfe et all Il99l [Storrie-Lombardi fc Wolfi ' 
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Noterdaeme et al.1 I2009T ). At lower redshifts, space- 
based spectra are necessary to measure the neutral 
hydrogen column densities in the UV. As DLAs (and 
QSOs that arc bright enough in the far ultraviolet 
(FUV) to be accessible to previous UV spectrographs) 
are relatively rare, the number of such systems currently 
known at z < 1 is small compared to available samples 
at high redshift. 

The Cosmic Orig ins Spectrograph (COS ) is a new in- 
strument package (jFroning fc Greenl I2009T ) installed on 
the Hubble Space Telescope that enables us to study 
these absorbers at low redshift with unprecedented ef- 
ficiency. As z < 0.5 spans ~ 40 percent of the age of the 
universe, the low redshift absorbers accessible with COS 
arc crucial for understanding cosmic chemical evolution 
and the cosmological gas mass density and for linking 
their properties to their higher redshift counterparts. 

In this paper, we report on the first DLA systems ob- 
served with COS. We present a detailed analysis of one 
such system, a DLA in the line of sight to the QSO SDSS 
J1009+0713, to illustrate the scientific potential of the 
observations. The structure of the paper is as follows. 
In §2 we describe in general the HST program and data 
reduction methods we have used for the COS spectra. 



In § 3 we describe the observations of the field of SDSS 
J1009+0713, including the ground based Keck/HIRES 
spectra of the QSO, the COS UV spectra, and imag- 
ing of the field with the Wide Field Camera 3 (WFC3). 
In § 4 we derive chemical abundances for this system, 
discuss the physical state of the g determined by 
the CI and C II* lines, and discuss the properties of the 
galaxies in the field as seen in the WFC3 images. In § 
5 we combine the measurements of N(H I) from other 
absorbers in this program to measure the cosmological 
mass density of neutral gas, fin i, at z < 0.35 in a blind 
survey with Az ~ 12. Conclusions are summarized in 
§ 6. Throughout this paper, we adopt a cosmological 
model with H m = 0.30, fi A = 0.70, and H o =70 km s" 1 
Mpc^ 1 . 

2. PROGRAM SUMMARY AND OBSERVATIONS 

The targets presented here were observed as part of 
program GO 11 598, a program focused on studying mul- 
tiphase baryons in the halos of L > L* galaxies at 15- 
150 kpc impact parameters. Target QSOs for this pro- 
gram were selected based on sufficient FUV flux and 
the presence o f a galaxy seen i n the Sloan Digital Sky 
burvey fSDSS. lYork et al.ll2000h with impact parameter 
p < 150 kpc and a spectroscopic or photometric redshift 
0.15 < z < 0.35. Sightlines with strong Mg II systems 
seen in the SDSS spectra were avoided due to the possi- 
bility of a Lyman-limit system (LLS) at z <0.35 which 
would block the FUV flux. In total, 39 sightlines will ul- 
timately be observed in this program. Followup ground- 
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Figure 2. The fit to the DLA line in the spectrum of SDSS J 1009+0713. The solid green line represents the adopted continuum and 
the red lines are the best fit model with log N=20.68 bracketed by profiles with column densities lcr=0.10 dex above and below the best 
fit value. The la flux uncertainty is shown in cyan below the data. 

ered in this program are shown in Figures Q] and [2] The 
Lyman alpha lines from the systems were fitted using 
a single component with the continuum and profile fit 
simultaneously. The redshifts of the systems were deter- 
mined from the metal lines, which with narrow features 
provide for a more accurate redshift determination. In 
total, we find three DLA systems and four sub-DLA sys- 
tems in these COS data. Absorption redshifts and neu- 
tral hydrogen column densities for the systems are given 
in Table [TJ Errors on Nh i given in Table Q] have been 
estimated by eye and are dominated by the continuum 
placement uncertainty. 

3. THE DLA IN SDSS J1009+0713 

In order to illustrate the capabilities of COS for fu- 
ture observations of DLA systems, we present a detailed 
analysis of one system observed in the sample here. Sim- 
ilar analyses for the complete DLA and sub-DLA sample 
will be given in a forthcoming paper. SDSS JI009+07I3 
(a=I0:09:02.0, <5=+07:13:43.8, z em = 0.456), is an op- 



based spectra with the Low Resolution Imaging Spec- 
trometer (LRIS) on Keck will provide for redshifts, star 
formation rates and metallicities of the target galaxies as 
well as galaxies in the field. 

These COS data were processed with CALCOS ver- 
sion 2.12 during retrieval from the archive. The indi- 
vidual reduced and wavelength calibrated xld files were 
then coadded using custom software. In short, the gross 
source and background counts in each wavelength bin are 
summed after the bins are aligned in wavelength space by 
using a cross correlation analysis to allow for small shifts 
due to inaccuracies in the wavelength calibration pro- 
vided by CALCOS. Simultaneously, an effective integra- 
tion time at each wavelength bin is created which the to- 
tal counts are ultimately divided by, returning the source 
count rate in photon s . During coaddition, the spec- 
tra are background subtracted and also corrected for the 
fixed pattern noise caused by grid wires above the micro- 
channel plate detectors using a reference file provided by 
the COS instrument team. COS has extremely low back- 
grounds, and consequently, in some cases there are very 
few counts in the cores of strong absorption lines. We 
determine uncertainties directly from th e accumulated 
counts, using Poisson statistics (see, e.g.. lGehrelslll986l ) 
when the count numbers are low. 

The spectra were binned by 3 pixels, as the raw COS 
data are oversampled with a <~ 6 pixel wide resolution 
element. All subsequent measurements and analysis were 
performed on the binned spectra. The binned spectra 
have a resolution of ~15 km s _1 per resolution element. 
The S/N of these COS spectra range from 7—15 per 
resolution element. 

The DLA and sub-DLA systems that we have discov- 



tically bright (m g =17.08, M fl 



-24.7) QSO observed 



m the SDSS. This QSO was chosen for the survey of 
galaxy halos due to the presence of an L ~ L* galaxy 
at an impact parameter of 60 kpc at redshift z = 0.228. 
The DLA described here was screndipitously discovered 
in the COS spectrum of the QSO. No absorption lines 
from this system are present in the SDSS spectrum of 
the QSO, as the Mg II AA 2798, 2802 lines which are 
typically strong enough to be seen in the SDSS spectra 
fall at 3115-3125 A and are below the wavelength range 
that the SDSS spectra cover. In this section, we present 
a detailed analysis of this system, including the COS and 
Keck High Resolution Spectrometer (HIRES) spectra as 
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well as followup imaging with the Wide Field Camera 3 
(WFC3). 

3.1. COS Spectra 

These COS data were acquired on 29-30 March 2010 
in three orbits as Visit 13 in program G011598. Two 
exposures were obtained with the G130M grating at cen- 
tral wavelengths of 1291 and 1309 A with exposure times 
of 1497 and 2191 s. Similarly, two exposures were taken 
with the G160M grating at central wavelengths of 1577 
and 1600 A with exposure times of 2002 and 2007 s. 

The DLA line of this system was fit with a single com- 
ponent centered at 2^=0.1140 with a best solution of 
log N(H I)=20.68±0.10. The fit to the DLA line is shown 
in Figure [2] Profiles with column densities ltr above and 
below the best fit value are also shown in Figure [2] 

Absorption lines from multiple species are shown in 
the left panel of Figure |()] We detect lines from N I AA 
1134, 1199, 1200, 1201; O I A 1302; Si II AA 1190, 1304, 
1526; S II AA 1250, 1253, 1259; Fe II AA 1142, 1143, 
1144; Ni II A 1317. A Lyman limit system (LLS) at 
z n .h.<i =0-3556 blocks the cont inuum blueward of ~1240 A 
(see iTumlinson et al.l ([2010) for an analysis of the LLS) 
or ~ 1115 A in the rest frame of the DLA. 

There are only two H 2 lines covered in the COS spec- 
trum due to the presence of the LLS. The Lyman band 
(X 1 !]^ — B 1 ^) transitions from the (1-0) band which 

have rest frame wavelengths A < 1109A are all blocked 
by the LLS. Only the H 2 AA 1120, 1125 lines resulting 
from the (0-0) transitions are covered in the COS data. 
These lines are however blended with strong high Ly- 
man series absorption lines from the LLS; consequently 
no direct indicator of molecular gas exists for this DLA. 

The COS spectra cover the lines of C IV AA 1548,1550 
as well as the Si IV AA 1393, 1402 transitions. The C IV 
lines lie at the edge of the detector where the sensitivity is 
low and only the C IV A 1548 line is minimally detected 
at S/N~3. The Si IV A 1393 line is blended with the 
Lyman-a absorption line from a z a ;, s =0.2774 system, and 
there is no detection of the Si IV A 1402 with S/N~5 in 
the region. 

3.2. Keck-HIRES Spectra 

Optical spectra of SDSS J1009+0713 wer e obtained 
with the Keck-II/HIRESb spectrograph (jVogt et al.1 
[1991 on the night of 26 March 2010. The S/N of the 
data ranges from <5 per 1.8 km s^ 1 pixel at 3200 A to 
^20 at 4500 A. These data provide coverage from 2750- 
4000 A in the rest frame of the z a&s =0.1140 DLA, and 
hence include the lines of Mg II AA 2796, 2803, Mg I A 
2852, Ti II AA 3242, 3384, and Ca II AA 3934, 3969. The 
Mg I A 2852 line is partially blended with the Fe II A 2344 
line from the z Q b,j=0.3558 LLS, as can be seen in Figure 
[H In addition, with higher resolution (R=44,000) than 
the COS spectra (R~20,000), these data provide valuable 
constraints on the kinematics and doppler parameters of 
the lines that are blended at the resolution of COS. 

These data were processed with the Keck/HIRES re- 
duction code hiresredux in IDlfl. The right panel of 
Figure [6] shows the HIRES spectra of several near-UV 

7 Available at http://www.ucolick.org/~xavier/HIRedux/ 



Table 1 

Properties of the z a ^ s < 0.35 DLAs and sub-DLAs in this survey. 



QSO 




z abs 


log N(H I) 








cm z 



SDSS J092554.70+400414.1 0.471 0.2477 19.55±0.15 

SDSS J092837.98+602521.0 0.295 0.1538 19.35+0.15 

SDSS J100902.06+071343.8 0.456 0.1140 20.68+0.10 

SDSS J143511. 53+360437.2 0.429 0.2026 19.80±0.10 

SDSS J155304.92+354828.6 0.722 0.0830 19.55±0.15 

SDSS J161649.42+415416.3 0.440 0.3211 20.60±0.20 

SDSS J161916.54+334238.4 0.471 0.0963 20.55+0.10 



Table 2 

AB magnitudes and position offsets for the galaxies near SDSS 
J1009+0713. Identifiers are based on the position angle and 
distance from the target QSO. 



ID 


Aa 


A<5 




m390W 


z gal 




// 


n 


AB 


AB 




170.9 


+1.58 


-9.27 


21.33+0.08 


22.16+0.12 


0.35569 


86.4 


+2.87 


+0.16 


23.21+0.21 


24.10+0.20 


0.35455 


80.5 


+4.00 


+0.81 


22.50+0.12 


23.09+0.18 


0.87899 


92.7 


+6.26 


-0.44 


23.67+0.37 


24.45+0.36 


1.283? 



(NUV) and optical absorption lines that were observed 
in this DLA. 

The profile was fit with only two components having 
doppler widths of b e ff — 10.6 and 12.9 km s _1 and cen- 
troid velocities of —6.6 and +21.5 km s _1 respectively. 
The majority of the low ionization gas appears to be in 
the v = —6.6 km s _1 component, as this component 
is strongest in all of the non-saturated absorption lines 
detected here. 

3.3. WFC3 Imaging and Keck LRIS Spectra 

The SDSS images show no obvious candidates for the 
z a &s=0.1140 DLA. In order to determine the nature of 
the host galaxy of the DLA system, we acquired imaging 
in the F390W and F625W filters of this field using the 
Wide Field Camera 3 (WFC3) onboard the HST. Six 
exposures were taken in both filters in a dither pattern 
to eliminate the effects of bad pixels and cosmic rays. 
Total exposure time was 2370 and 2256 seconds in the 
F390W and F625W filters respectively. 

Even with the superb image quality and resolution that 
WFC3 enables, galaxies at small impact parameter to the 
central bright QSO may only appear after subtracting 
the QSO with a suitably chosen Point Spread Function 
(PSF). The PSFs used here were acquired from the Tiny 
Tim online interface^ at the pixel coordinates appropri- 
ate for the position of the QSO in the images. This PSF 
was then interactively subtracted from each individual 
dithered image with idpE0 and the subtracted images 
were then drizzle combined with multidrizzle. The fi- 
nal drizzled images are shown in Figure [3] A color com- 
bined image was also created from the PSF subtracted 
images, and is shown in Figure H The F625W, F390W, 
and an average of the two were used in the red, blue, and 
green channels of the color image and th e images were 
scaled using the asinh scaling described in lLupton et al.l 
(|200l . 

8 Available at http : //www. stecf . org/ instrument s/TinyTim/tinytimweb 

9 Available at http://mips.as.arizona.edu/MIPS/IDP3/ 
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Figure 3. F625 and F390W band images of the field of SDSS 
J 1009+0713. The PSF subtracted images are shown in the bottom 
row. 
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Figure 4. Color combined PSF subtracted images of the field of 
J1009+0713 in the F625W and F390W filters. The red channel is 
the F625W band image, the blue channel is the F390W band image, 
and the green channel is an average of the two. 5 arcseconds is ~ 
10 kpc at z = 0.114. The arrow marks the location of the small 
feature visible after PSF subtraction. An H II region is marked 
with a box southwest of the center of the galaxy. North is up and 
east is to the left. 

Multiple galaxies are seen in these WFC3 images 
within 10" of the QSO, or - 20 kpc at z = 0.1140. 
A large spiral galaxy is clearly seen in Figures [3] and 
[3 Spectra of objects 80_5, 86_4, 92_7 and 170_9 as 
well as the H II region labelled in Figure 0] were ob- 
tained with Keck/LRl£3 • A 1" slit was used and aligned 
such that multiple galaxies and the H II region marked 
in Figure 2] were observed simultaneously. Emission 
lines from O III] AA 4960,5008 and H /3 A 4862 arc 

10 Labels for galaxies in the field are based on the position angle 
and distance from the QSO; i.e. galaxy 80_5 is at a position angle 
of 80 degrees east of north and 5 arcseconds from the center of the 
QSO. 



P625W 




Figure 5. Left - The WFC3 image of the QSO. Center - The 
model adopted from GALFIT. Right - The residual image. Some 
emission is left over in the point source in due to saturation effects. 
All frames have identical scaling. 

clearly seen in galaxies 86_4 and 170_9 placing them at 
z = 0.35458 and z = 0.35557 respectively. These two 
gala xies are likely associate d with a LLS at z a (, s =0.3556 
(see iTumlinson et al.l (|2010T ) for an analysis of the LLS). 

Two emission features from galaxy 80_5 are seen at 
7004 and 8158 A, which we interpret as O II] A 3727 and 
H 7 4341 at z = 0.87899. A single strong emission line 
at 8508 A is observed in the spectrum of galaxy 92_7. No 
prominent emission lines expected from a galaxy at a red- 
shift z = 0.1140 correspond to this observed wavelength. 
Assuming this emission line is from O II] AA 3727,3729 
places this galaxy at z = 1.283. 

Narrow emission lines from H j3 A 4862 and O III] AA 
4960,5008 from the H II region in the spiral galaxy la- 
belled in Figure U] give a redshift of z = 0.45699. Some 
flux from the nucleus of the QSO is observed in the spec- 
tra, however no emission lines from another redshift were 
seen in the spectra increasing our confidence in the de- 
termination of the redshift of the H II region. Combined 
with the close alignment of the nucleus of the QSO and 
disk of the galaxy, we interpret the spiral galaxy to be 
the host of the QSO itself. 

Emission lines from intervening galaxies are some- 
times seen superi mposed on the continuum of the QSO 
in SDSS spectra (iNoterdaeme. Srianand fc Mohanll20T0l : 
iBorthakur et al.ll2010f L although most QSOs do not show 
such emission features in their spectrum even when 
strong absorption lines are seen. No narrow emission 
lines are seen in the SDSS spectrum of this QSO at the 
absorption redshift of the DLA or any other redshift. 

4. RESULTS 
4.1. Imaging 

We performed aperture photometry on the galaxies in 
Figure|4]with APPHOT in I RAF. AB magnitudes are given 
in Table [2] Magnitudes h ave been corrected for Galac- 
tic du st extinction^ from iSchlegel. Finkbeiner &: Davisl 
(|1998h . As galaxies 170_9, 80_5, and 92_7 are spectro- 
scopically confirmed to be at significantly higher red- 
shift than the absorption redshift of the DLA, an obvious 
question arises; where is the host galaxy of t he DLA? 

We have used the galfit software package ([Peng et al.l 
|2010| ) to fit a multicomponent profile to the F625W im- 
age of the QSO. A PSF for the nucleus, Sersic pro- 
file for the bulge, and exponential disk for the spiral 
structure was used in the fit. The best fit solution 
from GALFIT using these components gives mg = 18.68, 
ni£> = 18.71, and mpf = 17.99 where tub, and rri/v 
are the magnitudes of the bulge, disk and nucleus respec- 
tively. The combined magnitude of the bulge and disk 

11 The online version of the calculator is available at 
http : //nedwww . ipac . caltech . edu/f orms/calculator .html 
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is mtot — 17.91, comparable to that of the nucleus. We 
show the results of the GALFIT decomposition in Figure 

El 

The host galaxies of QSOs are known to be luminous, 
with L>L*, and come from a variety of morphological 
types across the entire Hubble sequence dBahcall et al.l 
119971: iHamilton. Casertano fc Turnshekl 12002ft . Givena 
measured apparent magnitude m\, the absolute magni- 
tude M\ is given by 



AI\ = m\ - 51og 



(— 

V10 pc 



K A -A> 



(1) 



where Dl is the luminosity distance, K\ is the k- 
correction and A\ is the extinction due to the Galaxy. 
Based on the morphology of the galaxy seen in Fig- 
ure |4j we assume a Sb type spectrum and find a K 
correction of Ka=0.25 magnitudes in the F625W filter 
at z em = 0.456 and apply an extinction of A\=0.Q5 
magnitudes giving Mi?625W = —24.4. The median lu- 
minos ity of radio loud QSO host galaxies in the sam- 
ple of IHamilton. Casertano fc Turnshekl (|2002t ) is M v = 
—23.5, although galaxies classified as having spiral mor- 
phologies were typically more luminous than elliptical 
counterparts. 

A small feature at small impact parameter appears 
more clearly in the PSF subtracted images shown in Fig- 
ure [3] and is marked with an arrow in Figure 2] The 
angular separation from the nucleus of the QSO is 0.7 
arcseconds, or ~1.5 kpc at z = 0.114. It is possible that 
this is the host galaxy of the DLA, or perhaps an H II re- 
gion in the QSO host galaxy. For example, if either of the 
galaxies labelled 80_5 or 86_4 in Figure [4] were placed at 
small impact parameter to the QSO, it would be difficult 
to disentangle the two even with the resolution of WFC3. 
The WFC3 grism is likely the best (if not only) way to 
definitively determine the redshift the spiral galaxy and 
object at small impact parameter seen in Figure |U 

This QSO does show radio continuum emission, 
with a central source and two lobes detected in 
the FIRST survey with peak em ission of ~ 6 
mJy (jBecker. White fc Helfancj H991 . Radio loud 
QSOs are almost all found in elliptical galax- 
ies; out of the 26 radio loud QSO s observed by 
IHamilton. Casertano fc Turn shck (2002) , only 4 had spi- 
ral type morphologies (16 percent). 

4.2. Column Densities and Abundances 

Logarithmic abundances (defined as [X/H]=log X - log 
H -[X/H] Q ) for the DLA system as determined by the 
profile fits to the COS and HIRES data are presented 
here. Sola r system r eferen ce abundances ([X/H] Q ) are 
taken from lLoddersI (|2003[ ). Rest frame wav elengths as 
well as / values are taken from lMortonl (|2003l ) for all lines 
aside from the C I li nes, for which we use th e revised / 
values determined in Uenkins fc Trippl (|2001[ ) and Ni II 
AA 1 317, 1370 lines for which we use the revised / values 
from Uenkins fc Trippl (p006h . 

The COS spectra were fit using the component struc- 
ture determined with the higher resolution ground based 
Keck spectra. Both the dopplcr parameters and velocity 
centroids were held fixed in the fits, and multiple tran- 
sitions were fit simultaneously where possible (i.e. the 
S II AA 1250, 1253, 1259 lines were all fit simultane- 



ously). The synthetic Voigt profiles were convolved with 
the COS line spread functions (LSF)G3 as determined by 
the C OS science team at the nearest tabulated wave- 
length (jGhavamian et al.ll2009t ). The Keck/HIRES spec- 
tra were convolved with a gaussian line spread function 
with FWHM=2.5 pixels. Profile fitting was performe d 
using the profile fit code of lFitzpatrick fc Spitzerl (|1994| ) . 

In the case of a particular resonance line not being 
detected in these data, we have used the following pro- 
cedure to determine a limiting column density for the 
line. 1) First, we again use the doppler parameters and 
velocity centroids as determined in the higher resolution 
data. 3) We take the column densities for one of the 
better constrained fits to determine the ratio of the col- 
umn densities in each of the components. 3) Synthetic 
Voigt profiles are then produced and multiplied by the 
normalized flux and an equivalent width is determined. 
The integration and statistics are performed over a re- 
gion from the lowest to highest velocity component, with 
an additional 30 km s _1 (or ~2 resolution elements) on 
each end to ensure that all of the absorption is taken 
into account. 4) The significance of the line is deter- 
mined by calculating EW/ a tot where a tot = \J°~1. + °"p 

and a c is the uncertainty in the EW due to the con- 
tinuum placement and a p is the uncertainty This again 
perhaps indicates that the spiral seen in Figure |3] is not 
the host of the DLA galaxy, due to photon noise. 5) The 
column densities are then increased until the line has a 
3cr significance. 

Strong transitions such as O I A 1302, C II A 1334 and 
Mg II AA 2793, 2803 are almost universally saturated, 
and only lower limits to the column densities can be de- 
termined by fits to these lines. Lower limits in these cases 
were determined by creating synthetic Voigt profiles and 
increasing the column densities in the components un- 
til the profile reached the observed flux at the line core. 
We have determined a permissible range of 15.35<N(Mg 
II)<16.0 which is equivalent to -0.9 <[Mg/H]< -0.2 by 
measuring an upper limit on the weak transition of Mg 
II A 1239. 

As the component structure in the more highly ionized 
gas traced by the C IV A A 1548, 1550 and Si IV A A 1393, 
1402 lines is often different than that of the low ions 
and the S/N in the region of the C IV lines is too low 
to discern the component structure, we have measured 
the column density of C IV by direct integration of the 
line using the apparent optical depth (AOD) method. 
We determine log N(C IV)=13.86±0.20. We also place 
an upper limit on the column density of Si IV based on 
the non detection of the Si IV A 1402 line of log N(Si 
IV)<13.3. 

The depletion in this system is mild with 
S/Fe]=+0.24 ± 0.22, [S/Ti]=+0.28 ± 0.15 and 
L S/Ni]=+0.35 ± 0.22, even less than what i s see n 
in the halo of the Milky Way (|Weltv et al.l 119971) . 
Depletion levels in the ISM of the Milky Way are seen 
to vary over two orders of magnitude, from [Zn/Fc]> +2 
in cold clou ds to [Zn/Fe]~ +0 - 5 in t he halo of th e 
Milky Way (ISavage fc Sembachl [19961 : l.lenkinsl 12009ft . 
We note however that similar mild depletion was seen 

12 Tabulated LSFs for COS are available at 
http : //www. stsci . edu/hst /cos/per f ormance/spectral_resolut ion/ 
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Figure 6. Velocity plots for several lines detected in the COS and HIRES spectra towards SDSS J1009+0713 at z a6s =0.1140. The 
theoretical best fit profile is shown in red, and the la error in the COS spectra is shown below the data in cyan. The dashed vertical lines 
denote the positions of the centroids of the components used in the fits. Unrelated interloping features from other systems arc marked with 
magenta circles. Profile fits to nearby members in the multiplcts arc shown in the N I A 1134 and Fe II A 1142 lines. 



at very small impa ct parameter towards HS1543+593 
(jBowen et al.l l2005f ). a case in which the connection 
between the absorber and host galaxy is unambiguous. 

We have used the model of th e deplet i on pa tterns in 
the ISM of the Milky Way from iJenkinsI (|2009t) to ana- 
lyze the depletion in this DLA system. In this model a 
single value (F„) indicates the overall level of depletion 
in a sightline, with larger values of F* indicating higher 
levels of dust depletion. We find that the line of sight de- 
pletion factor F* = -0.48 ± 0.24 for this DLA, implying 
an intrinsic metallicity of [Z/H]=— 0.77±0.44, consistent 
with the metallicity determined from the S II lines. Such 
negative values of F* are typically seen along sightlincs 
with ex tremely low av erage densities in the halo of our 
Galaxy (|Jenkinsll2009f) implying that the sightline passes 
through the outer edge of a dwarf galaxy. 

Kincmatically, the component structure of this DLA is 
extremely simple (see Figure [6]) . Only two components 



Table 3 

Total column densities and abundances for the dominant ions in 
the z a (, s =0.1140 DLA system. 



Element 


log N(X) 


[X/H] 


H 


20.68±0.10 




C 


>15.90 


> -1.2 


N 


15.11±0.10 


-1.40 ±0.14 


O 


>16.00 


> -1.40 


Mg 


>15.35, <16.0 


> -0.9, < -0.2 


Si 


>15.00 


> -1.2 


1' 


<13.80 


< -0.34 


s 


15.25±0.12 


-0.62 ±0.16 


Ti 


12.70±0.03 


-0.90 ±0.10 


Fe 


15.29±0.17 


-0.86 ±0.20 


Ni 


13.93±0.18 


-0.97 ±0.20 



are required to fit the profiles of the absorption lines. 
The kinematic width Avqo (defined as the region where 90 



8 



Mciring et al. 



❖ 
o 




q This Work 
^ Pettini (3008) j 
a Van Zee (1998)- 
_ Van Zee (2006) : 



Figure 7. [N/a] vs. [a / Hi fo r this system (red point), high 
z DLAs from Pettini et al. (2008), H II regions in nearby dwarf 
galaxies from Van Zee & Havncs (2006), and H II regions in nearby 
massive spirals from Van Zee et al. ( 1998). The dotted lines denote 
the solar abundance. 

percent of the apparent optical depth in the line occurs) 
has often been used as an indi cator of the kinematics of 
QSO absorption lin e systems (|Prochaska fe Wolfd 119971 : 
iLedoux et al.l [20061) . We find a kinematic width for this 
system of Ai>go=52 km s _1 based on the Ti II A 3384 
line. 

4.3. Nitrogen 

The nucleosynthetic origins of N seem to include both 
primary (production of N from newly synthesized C via 
the CNO cycle) and secondary production (N produced 
again in the CNO cycle, but from previously created C 
and O from previou s star format ion) processes in inter- 
mediate mass stars (|Marigoll2001[ ). Of the ~30 determi- 
nations of N abundances in DLA systems, all are at high 
z where the N I lines are redshifted into wavelengths 
accessible with ground based telescopes ([Pettini et al.l 
120081 : iHenrv fe Proc haska 2007). Some constraints have 
been placed on nitrogen in low-z sub-DLAs where, like 
the high-z DLAs , Nitro gen is found to be undcrabundant 
(e.g. lTripp et~aLl IpOOll )). 

Abundances from N can be readily obtained from opti- 
cal spectra of H II regions in nearby star forming galaxies 
of a range of luminosities. Th e low luminosity galaxies 
from I Van Zee fe Ha vncs (2006|) show a flat plateau in the 
[N / O] ratio with respect t o metallicity, while th e higher 
metallicity H II regions in iVan Zee et al.l (|1998|) show a 
sharp increase in [N/O] with increasing metallicity, in- 
dicative of secondary production of N. 

Observations of N in extremely low metallcity DLAs 
are important for the study of N pr oduction and galac- 
tic chemical evolutio n in general (|Pettini et al.l 120081 : 
iCenturion et al.l |2003f). It is still yet to be seen if the 
[N/O] ratio in DLAs displays the same upturn as metal- 
licity increases as most DLAs are metal poor and at high 
z the N I lines are often blended with the Lyman-a for- 
est. The DLA presented here is the highest metallicity 
DLA with both N and a abundance measurements. 

In Figure [7] we show [N/a] vs [a /HI for H II re - 
gions in nearby galaxies from IVan Zee et al.l (|1998t ): 
I Van Zee fc Havnesl (120061 ). as well as the DLA sample 
from lPettini et al.l (|2008l ) and this DLA. The DLA from 
this work is shown with the red circular point, and lies 
below the "knee" of the diagram where N production 
transitions from primary to secondary pr oduction. The 
mean [N/a] ratio from the samples of IVan Zee et al.l 
(fl998l) and (|Van Zee fc Havnesl 12001 with [a/H]<-0.3 



is ([N/a]) = —0.55 ± 0.15, 3c x higher than seen i n this 
DLA of [N/S]=-0.78 ± 0.16. IPettini et al.l (120081) note 
this discrepancy in the sample of DLAs shown in Figure 
[7] and propose that the lower values seen in high z DLAs 
could be due to the delayed enrichment of N with respect 
to O as the delay of -~250 Myr is a greater fraction of the 
time available for star formation at high redshift than in 
local galaxies, so at high z this effect should be more pro- 
nounced. However, in this low-z DLA we see the same 
decrement in [N/a] compared to the values from local H 
II regions. It would be interesting to see if other high 
metallcity DLAs are also below the plateau, perhaps in- 
dicating a systematic offset in abundances determined 
via absorption and emission line diagnostics. 

4.4. C I and C 77* 

Both C I and C II* are important diagnostic resonance 
lines in the study of the ISM. Emission from the 158 fim 
emission line of singly ionized carbon in the ISM of the 
Milky Way is the dominant coolant, and is expected to 
also be so in high-z DLAs as well. This line originates 
from the transition between the 2 7-3/2 and 2 7\/2 levels 
in the ground state of C + . In the cold neutral medium 
(CNM, T~ 10 2 K) the excited state is populated largely 
by collisions with hydrogen atoms, where in the more dif- 
fuse warm neutral medium (WNM, T~ 10 4 K) both col- 
lisions with free electrons a s well as photon pumping can 
popu l ate the excited sta te (jWolfe. Prochaska fc Gawiserl 
[2001 iBowen et al.ll2005l) . 

In the WNM the electron density can be determined 
from observations of the C II and C II* lines as 



n(e) = 18.3. 



P 0.091/T 4 



N(cir) 



N(CII) 



cm 



(2) 



where T 4 is T/10 4 K (jLehner. Wakker fc Savagel 12001 . 
In the ISM of the Galaxy, the CNM is characterized 
by having large depletions of the refractory elements 
such as Fe and Ti relative to the volatile elements of 
S and Zn which do not condense onto dust grains. 
As the depletions levels in this DLA are modest, with 
[Ti/S]=— 0.37±0.12, we assume that the sightline is pass- 
ing through mainly a WNM. Depletion levels seen in the 
CNM of the Milky Way are typical ly ten times higher 
than what is obser ved in this DLA ([Savage fc Sembachl 
[l99i[Jenkmi[200l). 

The C II* A 1335 line was not detected in these 
data, with S/N^ll per resolution element in the region 
around the line. We have used the procedure discussed 
above to determine a limiting column density of log N(C 
II*)<13.7. We show a synthetic Voigt profile with this 
total column density overplotted on the data in Figure 

m 

As was mentioned above the C II A 1334 is strongly 
saturated and only a limiting column density of log N(C 
II)>15.90 can be determined from this transition. How- 
ever, in sightlines with low depletion such as this Carbon 
is mildly depleted and we can use the S abundance to es- 
timate a total C II column density assuming [C/S]~0.0. 
The above analysis yields log N(C II)=16.45. Using the 
above values in Equation [2] with a temperature of 6000 
K characteristic of the WNM the electron density is de- 
termined to be n(c)<0.08 cm -3 . 

The cooling rate per H atom due to C II] A 158 fim 
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Figure 8. Synthetic Voigt profiles of the CII* AA 1335.6, 1335.7 
lines overplotted on the COS data. We estimate the limiting col- 
umn density of CII* to be log N(C II*)~ 13.70. 




Figure 9. Synthetic Voigt profiles of the CIA 1277 line over- 
plotted on the COS data. We estimate the limiting column density 
of CII* to be log N(C II*)~ 13.6. 



emission can be estimated by 
N(C II*)- 



2.89 x 10 



-2(1 



N(H I) 



erg s 1 (H atom) 1 (3) 



(|Lehner, Wakker fe Savage! 12004ft . Using the val- 
ues above, we determine log I c < —26.5. This 
is typical of wha t has been observed in h igher- 
z DLA systems (IWolfe Prochaska fc Gawiseil 120031 : 



ILehner. Wakker fc Savage!l2004ft 



C I is a tracer of cold gas in t he ISM, and is usually seen 
in DLA absorbers bearing H2 l)Srianand et al.ll2005ft . As 
such, C I absorption lines are indicative of cold, dense 
material in DLAs. The CIA 1277 line provides the 
strongest constraint on the C° column density in the COS 
spectrum as the S/N in the region around the stronger 
transition of C I A 1560 was much lower. This line was 
not detected in the spectrum and we have determined 
a limiting column density of log N(C I)<13.6 usin g the 
prescr iption above. In the two phase ISM model of iLisztl 
pOOl . the ratio of N(C II)/N(C I) is indicative of the 
gas density. The observed ratio of N(C II)/N(C I)>700 
in this system implies that 11(H) < 10 cm -3 . 

5. THE COSMOLOGICAL NEUTRAL GAS MASS DENSITY 

The sample of QSO sightlines from this HST program 
provides an opportunity to determine the incidence of 
DLA systems at low z, where there is a paucity of data 
due to the necessity of space-based observations. Ideally, 
a totally "blind" survey in which there is no prior infor- 
mation about the sightlines prior to observation provides 



for the least amount of systematic biases in the determi- 
nation of cosmological gas mass density. 

In order to eliminate as much observational bias as pos- 
sible, we have defined the redshift interval over which to 
search for DLAs and sub-DLAs with the following crite- 
ria. As these sightlines were selected due the presence 
of a galaxy within impact parameters of p < 150 kpc, 
we have excluded ±1000 km s _1 around the redshift 
of the target galaxy in the search path to avoid biasing 
the sample artificially high due to the possible clustering 
of galaxies. At z > 0.34 the Mg II AA 2796, 2803 lines 
are redshifted into the wavelength coverage of the SDSS. 
This program intentionally excluded sightlines with Mg 
II absorbers at z > 0.34 due to the possibility of a LLS 
which would block the FUV flux, which in effect biases 
the sample against DLA detection at z > 0.34. To avoid 
this bias, we focus only on z < 0.34. In the cases where 
the emission redshift of the QSO is less than 0.34, we 



have excluded a 3000 
emission redshift. 

The sub-DLAs 
J0928+6025, and 
eluded from the 
shifts are within 
the target galaxy, 



km s region blueward of the 



towards SDSS J0925+4004, SDSS 
SDSS J1435+3604 have been ex- 
calculation as the absorption red- 
1000 km s" 1 of the redshift of 
leaving 3 DLAs and 1 sub-DLA 



in the final sample. The total path length for 
the 37 sightlines in this sample using the definition 
above is Az=11.94, giving d/V jdz = 



0.25lg;^ and 



d/V /dz=0.08Io 5g for DLAs and sub-DLAs resp ectively, 
and a ssuming Poisson statistical errors from iGehrelsl 
(|1986ft . The redshift density determined here is larger 
than the value determined in iRao. Turnshek fc Nestorl 
(|2006l) of d/V /dz=0.079±0.019, however the large sta- 
tistical errors of this work make this difference un- 
likely to be statistically signi ficant. From a large num- 
ber of archived HST spectra, iRibaudo. Lehner fc Howkl 
(|2010| ) find (W /dz=0.37±0.10 for all absorbers with N(H 
I)>17.0. The central value determined here for DLAs 
alone is ~2/3 of the total redshift density of all optically 
thick systems, again suggesting that the redshift density 
of DLAs is lower than what is determined here. 
The comoving gas mass density can be determined as 



fi H i = 



pm H H EA H : 
AA 



Cp c 



(4) 



where p = 1.3 is the mean molecular mass of the 
gas taking into account the contribution of He, p c is 
the critical density of the Universe, Ho is the Hubble 
constant, and mg is the mass of the hydrogen atom 
(| Wolfe. Gawiser fc Proc haska 20051). The absorption dis- 
tance dX is defined as 



dX 

dz 



(l + zf 



y/n m (i + zf + n k (i + zf + n A 



(5) 



and is summed over each sightline in the sample. The 
total path length for the 37 QSOs in this sample using 
the criteria above is AX=15.1. 

We find H i =l-*to% x 10" 3 and H 1 = 
10~ 5 for DLAs and sub-DLAs respectively. The neu- 
tral gas density in DLAs at z ~ 0.6 has been deter- 
mined by previous surveys to be Oh i = (9.7±3.6) x 10~ 4 
(|Rao fc Turnshek! [20001: IRao. Turnshek fc Nestorl 12006ft . 
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The value determined here agrees within the errors with 
these previous determinations. The low-z mass density 
of neutral gas from sub-DLA syste ms has not yet been 
thoroughly investigated. At z ~ 2.0. lPeroux et al.l f|2005l ) 
find H i =(1.9 ±1.3) x 10~ 4 . 

Due to the small number statistics, this sample is dom- 
inated by statistical errors. Upper and lower bounds on 
the gas density have been determined by using the up- 
per and lowe r ltr li mits on the number of systems, n , 
fromlGehrelsj (fl986h assuming Poisson statistics. We es- 
timate the la range of E N(H I) in Equation|4]as E N(H 
I)=?i ± (N(H I)} where (N(H I)) is the mean observed col- 
umn density in the sample. 

6. SUMMARY AND DISCUSSION 

In this work, we have reported on the first Damped 
Lyman-a systems observed with the Cosmic Origins 
Spectrograph. We summarize the results of this work as 
follows. 

1) We have determined the redshift density and gas den- 
sity of DLAs and sub-DLAs in this sample covering a 
total redshift path of Az = 11.9 at z < 0.34. We 
find <W/dz = 0.25±°;^ an d H i =1.4±J;| x 10~ 3 for 
DLAs, and for sub-DLAs dJV /dz=0.08±g;^ and H i 

=4.2±|;| x 10" 5 . The redshift density of DLAs deter- 
mined here is several times higher than previous measure- 
ments of this value at z ~ 0.6 of dA/"/dz=0.079±0.01 9 
(jRao fc Turnshekl [20001: iRao. Turnshek fc Nestor! I2006D . 
The large statistical errors of our sample make any dif- 
ferences unlikely to be statistically significant. No mea- 
sure of cUV / dz or fin i exists for sub-DLAs at z < 
1.5, however the value determined here is lower than 
dAf/dz= q.34±g;^ and Q H i = (1.9 ± 1.3) x 10" 4 seen at 
z ~ 2.0 (|Peroux et al.l 120051) . The neutral gas density 
in DLAs measured here is slightly higher than what has 
been found in previo us surveys at z ~ 0.6 which find Oh i 
= (9.7 ±3.6) x 10~ 4 (IRao. Turnshek fc Neston[2006h . but 
the allowed range is fully consistent with the previous 
measurements. Based on 2 1 cm emission map s of galax- 
ies in the nearby Universe, (jZwaan et alJ feOOS) find Oh i 
= (3.5 ± 0.8) x 10 -4 , which is well below the value deter- 
mined here. 

2) We have measured the abundances of several elements 
of one DLA in the sample and shown that the system has 
a sub-solar total abundance with [S/H]=— 0.62 ± 0.18. 
This metallicity is significantly higher than typical 
high redshift DLAs, which have average mctallicities 
([Zn/H]) ~ —1.25. The kinematical structure of the 
DLA is remarkably simple. Only two components are 
needed to fit the profiles and the velocity width is 
well below that of the average DLA population with 
Ai>go=52 km s _1 . Remarkably, the kinematics are 
similar to tho s e of th e extremely metal poor DLAs from 
iPettini et all (|200l) . even though the metallicity of 
this DLA is ^100 times greater than the DLAs in that 
sample. 

3) Followup imaging with WFC3 reveals several star 
forming and presumably gas rich galaxies within 10" of 
the QSO, as well as a large spiral galaxy at small impact 
parameter to the bright QSO. At first glance, any of 
these candidates appear to be promising candidates 



for the host galaxy of the DLA system. Followup 
spectroscopy of the galaxies reveal that none of them 
are at a similar redshift to the DLA. Spectra of an H II 
region in the spiral galaxy show strong narrow emission 
lines at the redshift of the QSO. This, combined with 
the close alignment of the nucleus and disk of the galaxy 
leads us to identify the spiral galaxy as the host of the 
QSO. The QSO host is luminous, with M F625 w = -24.4 
based on a decomposition of the WFC3 image with 
GALFIT. A small feature is seen ~0.7 arcseconds from 
the center of the QSO nucleus after PSF subtraction, 
adding another host galaxy candidate of the DLA. 
Even with this supporting data, the origin of the DLA 
host galaxy remains enigmatic and exemplifies some of 
the difficulties of discovering the host galaxies of such 
systems. 

4) Independent of the morphological properties of 
the DLA host galaxy, we conclude that the sight- 
line likely passes through the warm neutral medium 
of the host as the metals are very mildly depleted 
with [S/Fc]=+0.24 ± 0.22, [S/Ni]=+0.35 ± 0.22 and 
[S/Ti]=+0.28 ± 0.15 and both C II* and C I are 
not detected. The depletion levels and low cold 
gas co ntent are charact e ristic of higher z DLA sys- 
tems dProchaska et all l2002t iMeiring et all 120061: 
iNoterdaeme et al J 1200*91) . We find that the line of sight 
depletion factor F* = —0.11 ± 0.24, characteristic of 
sightlines through the halo of the Milky Way, reaffirming 
the conclusion that the sightlinc is passing through warm 
neutral gas. 

In a forthcoming paper we will report on the metallic- 
ities and depletions in the other DLAs and sub-DLAs in 
this sample, as well as the N(H I)-weighted mean metal- 
licity at z < 0.34. Few observations of S or Zn, un- 
depleted elements in the ISM and accurate metallicity 
indicators, exist at z < 0.5 which spans ~ 40 percent 
of the age of the universe. No observations of Zn or S 
in sub-DLA systems exist in this redshift range. At low 
rcdshifts, where there has been more time for the suc- 
cessive buildup of metals in the ISM of galaxies and the 
abundances of DLAs should be higher than at high red- 
shift, it will be possible to see if DLAs show an increase 
in the [N/a] ratio characteristic of secondary production 
of N. Our data show that over its lifetime, the Cosmic 
Origins Spectrograph will have the ability to measure the 
abundances and physical properties of QSO absorbers at 
z r~i with nearly the same precision as what can be done 
at z ~ 2 with large ground based observatories. 
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